Dipolar Stabilization of an Attractive Bose Gas in a One Dimensional Lattice 
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We experimentally show that dipolar interaction can stabilize otherwise unstable many-body 
systems like an attractive Bose gas. In a one dimensional lattice the repulsive dipolar on-site 
interaction balances negative scattering lengths up to —17 Bohr radii and stabilizes the ®^Cr Bose- 
Einstein condensate. For reduced lattice depths, the dipolar stabilization turns into destabilization. 
We probe the full cross-over between the two regimes and our results are in excellent agreement 
with theoretical calculations, which reveal significant dipolar inter-site interactions. 
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Recent years have seen an increasing interest in the 
design and realization of novel quantum phases. In this 
regard, cold quantum gases play a central role [l|, Q 
through the high level of experimental control that they 
allow. Internal as well as external degrees of freedom 
are controlled independently and one is able to tunc the 
two-body short-range interaction that drives the main 
properties of such cold ensembles almost at will Q . 

The toolbox for the design of possible novel quan- 
tum phases has been recently widened by a new ingre- 
dient: the dipolar interaction (DI). In strong contrast 
with the contact interaction, it has two assets, namely 
the anisotropy and the long-range character. In the last 
years, the DI in cold gases has been object of numer- 
ous experimental and theoretical investigations [1, 
Within this framework, quantum gases with strong dipo- 
lar interactions have been experimentally realized in sys- 
tems of chromium Bose-Einstein condensates (BECs) [6|- 
@] , ultra-cold heteronuclear molecules Q and BECs with 
Rydberg excitations [l^ , while weak dipolar effects have 
also been observed in alkali samples 11 - 3] ■ On the the- 
oretical side, unique self-organized structures have been 
predicted, such as vortex lattices of different symme- 
tries density modulated ground-states 15, 11], and 
supersolid phases in optical lattices 17- loj. 



When confined in a quasi-two-dimensional (2D) ge- 
ometry, the underlying mechanism responsible for self- 
organization in dipolar quantum gases is closely related 
to their roton-maxon excitation spectrum (20| . Further- 
more, in a stack of such quasi-2D systems, long-range 
inter-site interactions are responsible for an enhancement 
of dipolar effects [U, [23 and the amplification of the 
self-organized structures [23| . On the level of many-body 
physics, the dipolar atomic gases already fulfill the re- 
quirements of degeneracy and strong dipolar interactions. 
However, in these systems the structures described in 
Refs. 

are expected to form in a so far unreached 
range of parameters, where repulsive dipolar interaction 
counterbalances attractive contact interaction. Therefore 
investigating the stability of multi-site systems of highly 
oblate dipolar condensates is a crucial step towards the 



realization of novel quantum phases. 

Experimentally, a one dimensional (ID) optical lattice 
oriented along the polarization direction of the dipoles al- 
lows for the realization of the desired oblate BECs with 
mainly repulsive DI. One expects the clouds in the indi- 
vidual lattice sites to be stable even at negative scattering 
lengths due to a strong dipolar stabilization Q. But, as 
several lattice sites are populated, one also expects the 
mainly attractive inter-site interactions due to the long- 
range dipolar potential to destabilize the system. Be- 
cause of this interplay between on-site stabilization and 
inter-site destabilization, it is a priori unclear whether 
one can reach the regime of attractive contact interaction 
where different new quantum phases have been predicted 
to appear. 

In this work, we show that in a ID optical lattice the 
strong dipolar repulsive on-site interaction can, in fact, 
stabilize a BEC of ^^Cr atoms with attractive short-range 
interaction. We identify the stability threshold of the 
condensate for various lattice depths by measuring the 
BEC atom number when decreasing the s-wave scatter- 
ing length via magnetic Feshbach resonance. From small 
to large lattice depths we observe a continuous cross-over 
from a dipolar destabilized to a dipolar stabilized regime 
where we find a stable BEC even at negative scattering 
length. Indeed for deep lattices one can consider our sys- 
tem as a stack of mesoscopic quasi-2D BECs containing 
up to 2000 atoms in the central wells, strongly stabilized 
by the on-site dipolar interaction. Our measurements 
arc in good agreement with numerical mean-field calcula- 
tions, which include significant contribution of inter-site 
interactions mediated through the long-range dipolar po- 
tential. 

Our experimental procedure is as follows: we produce a 
^^Cr BEC in a crossed optical dipole trap (ODT, i^x.y.z = 
(440, 330, 290) Hz) at a magnetic field B ~ 600 G, where 
the scattering length is large and positive (a ~ 90 oq with 
ao the Bohr radius). Dipoles are aligned along the strong 
magnetic field in z-direction. By changing the magnetic 
field strength B in the vicinity of a Feshbach resonance 
(FR), we tune the s-wave scattering length a according 
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the determination of the stabihty threshold. 
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Figure 1: Experimental setup. The measurements are per- 
formed in a ID optical lattice (blue) with underlying crossed 
optical dipole trap (red) . The magnetic field used to reach the 
Feshbach resonance is produced by Helmholtz coils (black) 
and polarizes the dipoles along the lattice direction (z). For 
deep lattices we obtain a stack of oblate dipolar BECs as de- 
picted on the lower right. 
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where Obg — IOOoq is the background scattering length, 
A ~ 1.6 G the width and = 589.1 G the center of the 
FR. After reducing the scattering length to a = 60 Oq, 
we load the BEG into the ID optical lattice which is ori- 
ented along the polarization direction z. The lattice is 
produced by a A = 1064 nm fiber laser in a nearly back- 
reflected geometry (full crossing angle a = 10°) with a 
lattice spacing diat = 534 nm, as illustrated in Fig.[TJ The 
radial trapping frequencies {vx,y) are kept constant dur- 
ing the lattice ramp by adjusting the power in the ODT 
beam in the z-direction. We then decrease the scattering 
length in 6 ms to reach its final value, where we hold the 
atoms for thoid = 2 ms. Finally we switch off the optical 
trapping potential for a 6 ms timc-of-fiight (TOF) before 
taking an absorption image. 

To extract the BEG atom number iVeEC after TOF, 
the recorded 2D density distribution is integrated along 
the ^;-direction and we perform a ID bimodal fit. When 
the final scattering length is much larger than the crit- 
ical value flcrit, we typically measure A'bec — 15.000 
(-/Vbec — 10.000) in a shallow (deep) lattice, while 
the atom number before loading the lattice is always 
A^at — 20.000. Getting close to the instability point, we 
observe a fast decrease in the atom number as it is shown 
in Fig. [2] for two different values of the lattice depth. 
We finally extract the critical scattering length from an 
empirically chosen function as described in Ref. Al- 
though atom losses are enhanced in a deep lattice due to 
the larger mean trapping frequency, they do not affect 
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Figure 2: Atom number versus scattering length in shal- 
low and deep lattice. In a rather shallow lattice of depth 
U = (7.5 ± 0.8) -Br (open blue dots), the condensate becomes 
unstable at Ocrit = (6.5±1.9)ao, while in a deep lattice a.tU = 
(45 ± 5) -Er (filled red dots) we even find a stable BEC until 
a = ( — 13.2 ± 2.5) Oq. The solid lines are fits to the data using 

the arbitrarily chosen form A'bec = max |o, A^o • (a — flcrit)'' | , 

from which we extract the critical scattering length Ocrit (typ- 
ically f3 ~ 0.2). 

Figure [3] shows the stability diagram of a dipolar ^^Cr 
BEC in a ID optical lattice. The critical scattering length 
flcrit is measured for different lattice depths in the range 
[/ = - SOSr (recoil energy, -Er = H'^tt'^ /{2mdf^^), with 
m the atomic mass). We find positive a^it until U ~ 
10 and down to acrit = (— 17 ± 3)ao in a deep lattice 
potential. The datapoints are in very good agreement 
with numerical simulations based on the non-local non- 
linear Schrodinger equation 
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where g = 47rfi,^a(i?)/m. The potential t4xt(r) = 
f7sin^(7rz/diat) + mJ2i=x y z ^Tr^t'f ^1/2 results from the 
ID optical lattice and the 3D harmonic confinement given 
by the ODT. The DI potential is given by Vdd(r) = 
^^f^ '^"^lli^'^ (with f = r/r), where no is the vacuum per- 
meability and ^ the magnetic dipole moment (/i = 6^b 
for ^^Gr with /j,b the Bohr magneton). 

The DI term is calculated by using the convolution 
theorem and the Fourier transform of the DI potential, 
V;id(k) = 2i£|ii! 1^1 _ 3|^ X zpj [|. We evaluate the 

ground-state of the system using a full 3D simidation of 
Eq. ^ in imaginary time. The full 3D character of the 
calculations allows us to determine the critical scattering 



3 



length Ccrit for all lattice depths, ranging from U ~ to 
very deep lattices. Alternatively, we evaluate Eq. ([2]) in 
real time, and simulate the actual experimental sequence 
described above. The time-dependent results are basi- 
cally undistinguishable from those obtained from ground- 
state considerations in imaginary time. This implies that 
the experimental results recover the stability threshold of 
a dipolar BEC even for deep lattices, for which reaching 
experimentally the ground state becomes questionable. 
Furthermore, the comparison of simulations for different 
atom numbers (Fig. [3]) shows that the stability threshold 
depends only weakly on the initial BEC atom number 
A^at, in particular in the deep lattice regime. 
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Figure 3: Stability diagram of the dipolar condensate in tiie 
ID optical lattice. The critical scattering length acrit(ao) is 
plotted versus the lattice depth [/(_Er). The lines are results 
of the numerical simulations for different atom numbers. The 
full cross-over from a dipolar destabilized (flcrit > 0) to a 
dipolar stabilized (acrit < 0) regime is observed. At [/ « 
10 Eb. a purely dipolar interacting BEC (acrit = 0) can be 
stabilized. The dash dotted line (green) shows the simulated 
critical scattering length disregarding dipolar interaction. 

In order to understand the results shown in Fig. [31 
we recall that the stability of a dipolar BEC strongly 
depends on its geometry [7|. It is also convenient to 
introduce a length scale associated to the DI, add = 
TO/io^^/127r?i^ (~ 15 ao for ^^Cr). When the lattice is 
absent, the ODT determines the shape of the BEC and 
therefore its critical scattering length. In our case, the 
trapping frequencies arc chosen such that both trap as- 
pect ratios A^^^j^ = -jj^ are smaller than unity. In this 
essentially prolate trap, the dipolar interaction strongly 
destabilizes the condensate. The BEC becomes then un- 
stable at a positive scattering length acrit = (12 ± 2) oq, 
close to the value a = add, which is the expected 
critical scattering length in a prolate trap in the limit 



A^atOdd/oho ^ 1 [3 1 with aho the mean harmonic oscilla- 
tor length. 

On the contrary, for a deep lattice, the relevant aspect 
ratio is that of the single site trap, which becomes largely 
oblate. Hence, for a growing lattice depth, we observe a 
smooth cross-over from a dipolar destabilized (acrit > 0) 
to a dipolar stabilized (ccrit < 0) regime, characterized 
in Fig. [3] by a rapid decrease of acrit for depths J7 ^ 5 — 
10 -Br. For larger depths, at which inter-site hopping is 
negligible, the decrease of acrit becomes much slower and 
we find a minimum critical scattering length of acrit = 
(— 17 ± 3) ao at J7 = 60 Br. This is in strong contrast 
with the behavior of a purely contact interacting BEC 
(with same parameters otherwise), for which simulations 
(see Fig. |3]) show only weak dependence of acrit on the 
lattice depth, with acrit — Oao in the full range [j^ . 
In this respect, we show here that the DI stabilizes an 
attractive Bose gas. 

In the deep lattice region {U > 15 Br), a naive picture 
based on the stability criterion for a single oblate con- 
densate would suggest that acrit should rapidly approach 
— 2add- Experimentally, we observe a stable condensate 
until acrit — 17ao only. Indeed in a deep lattice, our 
system becomes a stack of quasi-2D pancake-like conden- 
sates with inter-site interactions mediated by the long- 
range DI. Due to the anisotropy of the DI, the inter-site 
interactions have a destabilizing character, and as a re- 
sult acrit is considerably higher than the expected — 2add- 
Note that this inter-site destabilization requires an in- 
homegeneous BEC, since two parallel infinite homoge- 
neous quasi-2D condensates would present zero-averaged 
inter-site DI El. 
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Figure 4: Inter-site coupling mediated by dipolar interactions: 
zoom on Fig. [3] in the deep lattice regime. Inter-site hopping 
is negligible on experimental time scales. Solid line: full 3D 
numerical simulation. Dashed line: simulation with the trun- 
cated dipolar potential (see text), for which inter-site cou- 
pling is not taken into account. The deviation of the simu- 
lation with the truncated potential to the experimental data 
indicates significant mean-field energy contribution from long- 
range inter-site interactions. 
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In the following, we further investigate the desta- 
bilizing role of the dipolar inter-site interactions. 
In this respect, strong insight is gained by in- 
troducing a truncated DI potential, Kid°'^('^) ~ 
Fdd(r) [e ((r • z) -t- diat/2) - e ((r • z) - d,at/2)], where 
0(^) is the Heaviside function. Such cut-off is im- 
plemented following a similar procedure as that of 
Ref. [29|. For deep lattices the cut-off in V^^'^iv) effec- 
tively amounts to remove inter-site DI, while still taking 
into account the short-range and long-range on-site inter- 
actions. In Fig. [5] we compare the stability threshold cal- 
culated with Vjd™('^) *o data and the full simulations 
in the deep lattice regime {U > 15 -Br). The calculations 
with the truncated dipolar potential show a substantial 
deviation to the experimental data (s^j: for instance, at 
U ~ 20 Ey{ the difference is Aocrit — 8ao, which is more 
than three times our standard deviation. This discrep- 
ancy, in addition with the very good agreement of the 
full simulations with the measurements, provides a clear 
signature of the presence of inter-site interactions. Note 
that these inter-site interactions reduce the stabilization 
effect of the on-site repulsive DI (e.g. for U = 20£'r, 
[acritl is reduced by close to 50%), indicating that they 
contribute significantly to the overall energy of the dipo- 
lar BEC in our experiments. 

In conclusion we have experimentally shown the sta- 
bilization of a polarized dipolar Bose gas with attractive 
contact interaction. We produced the required large trap 
aspect ratio by loading the ^^Cr BEC into a one dimen- 
sional optical lattice potential. The measurements are 
well described by numerical mean-field calculations us- 
ing the non-local non-linear Schrodinger equation. In the 
regime of negligible inter-site hopping, the data supports 
significant long-range inter-site interactions. 

The successful stabilization of a dipolar quantum gas 
with attractive short-range interactions paves the way to 
the investigation of the roton-maxon excitation spectrum 
expected for such systems 2^, 21, Slj and the related self- 
organized structures. Furthermore, the realization of a 
stack of stable mesoscopic ensembles interacting through 
long-range dipolar interactions, together with the pos- 
sibility to tailor multi-well potentials at will [13, Hjl , is 
particularly promising for the experimental realization of 
new quantum phases relying on the interplay between on- 
site and long-range inter-site interactions, as predicted in 
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